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Abstract 



The scaling limit of the higher level Bethe Ansatz (HLBA) equations for 
a macroscopically half-filled Hubbard chain is considered. These equations 
practically decouple into three disjoint sets which are again of the BA type, 
and correspond to the secular equations of three different kinds of dressed 
particles (one massive and two massless). The finite size corrections and the 
fine structure of the spectrum show that the massless sector corresponds to a 
conformal field with central charge c = 1 and Gaussian anomalous dimensions. 
The zero temperature free energy is also calculated and is found to be in 
perfect agreement with the results of a perturbative calculation in the SU(2) 
chiral Gross-Neveu (CGN) model. Some further arguments are presented 
supporting the identification of the model obtained as the relativistic limit of 
the half-filled Hubbard chain with the SU(2) CGN model. 

PACS numbers: 05.30.Fk, 05.50.+q, ll.10.Kk, 75.10.Lp 
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I. INTRODUCTION 



The one dimensional (ID) Hubbard model, being completely integrable |TJ and exactly 
solvable by Bethe Ansatz (BA) plays a central role in the theory of strongly correlated 
electron systems 0. It is described by the Hamiltonian 

N N / 1\ / 1\ 

# = E (i^v+iiA +^Ekt-« r^"2 ' (L1) 

j=i<t=u i =1 V 7 V 7 

where cf a resp. c J cr are the canonical electron creation resp. destruction operators at site 
j with spin a, and hj >a = c- a Cj a . We assume periodic boundary conditions, so the site 
j = N + 1 is equivalent to the site j = 1. The sign of the hopping, t, is positive while the 
interaction, U, for the considered attractive case is negative. The Hamiltonian ( |1.1| ) has an 
SU(2)xSU(2) symmetry, which is in fact reduced to SO(4) due to the selection rule that 
(for even N) only representations with integer total spin of SU(2)xSU(2) are allowed 
One of the SU(2) is associated to the spin the other we call isospin (which is associated to 
the charge). 

The diagonalization of the Hamiltonian ( |1.1| ) is reduced by the BA to solving a set of 
nonlinear equations 0. For a macroscopically large number of electrons the corresponding 
system consists of a macroscopically large number of equations, and to solve them means to 
reduce the original set to a smaller set which describes the excitations only. This reduced 
set is often referred to as the higher level Bethe Ansatz (HLBA) equations, and can be 
considered as the secular equations for the excitations. In this procedure also the dispersion 
of the excitations is found. The different solutions of the BA equations corresponding to 
different kinds of eigenstates of the Hamiltonian ( |1.1| ) are discussed in detail in the papers 
reprinted in Ref. |§. 

A macroscopically half-filled (HF) Hubbard chain (i.e. in which the number of particles 
does not significantly differ from the number of sites) has two types of excitations f7j|| which 
are connected with the two SU(2) symmetries of the system ||. One kind of the excitations 
has a gap in the energy spectrum, the other is gapless. The energy and the momentum of 
the system is given as the sum of the energies and the momenta of the individual excitations. 

In the present work we study the excitations (i.e. the physical or dressed particles) in 
the scaling limit of the HF Hubbard chain. To gain insight into the structure of the excited 
states, i.e. into the nature of the physical particles of the limiting theory, we derive their 
secular equations. The starting point of our analysis are the HLBA equations derived in 
Refs. and we follow the behaviour of these equations and the spectrum in the scaling 
limit. The scaling limit corresponds to a continuum limit together with fine tuning the 
interaction strength U to zero while keeping the mass of the physical particles described by 
the limiting HLBA eqs. fixed. In the continuum limit the lattice spacing, a, tends to zero 
keeping the length of the chain L = Na to be constant (so iV — > oo), while the particle 
number per site is fixed. In order to keep the mass scale (the energy gap) finite as a — > one 
has to adjust the hopping t together with the dimensionless coupling u = \U\/4t according 
to 

m o = — Vuexpl- — ), t(— l/2a) — > oo, u — > , (1.2) 



7T 2m 
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in agreement with the result found in Refs. PJT0||. As the spectrum of the Hamiltonian 
( |1 . ID in this limit contains both massive and massless excitations with relativistic dispersion 
relations, it is very tempting to identify the continuum model defined by the scaling limit 
with a (l+l)-dimensional relativistic field theory. 



It has been argued already a long time ago [II] that the scaling limit of the HF Hubbard 
chain can be identified with the SU(2) symmetric chiral Gross-Neveu []E| (CGN) model. 
The Lagrangian of the CGN model can be written as: 

C = ijfaf, - -g^a^) 2 - \g'{^) 2 , (1.3) 

where ip is a doublet of Dirac fermions, and a % are the Pauli matrices. This theory is asymp- 



totically free |L3| and it belongs to the very special class of field theories which admit an 
infinite number of conservation laws preventing particle production in scattering processes. 
Using the 'bootstrap' approach the .S'-matrix of the massive particles had been proposed 



11 16 1 . Furthermore the CGN theory is expected to essentially decouple into a massive 



and massless sector |T7[. The Hamiltonian of ( |1.3| ) has been diagonalized by the BA ]TB| , |T§ 
Thereby both the proposed mass spectrum and the bootstrap ^-matrix has been confirmed 



In Ref. |10[ both the spectrum and the S matrix of the massive particles computed in 
the scaling limit of the HF Hubbard chain has been found to agree with those of the SU(2) 
CGN model. The S'-matrix of the massless sector was found to correspond to the massless 
scattering theory proposed in Ref. |2]J associated to a level one SU(2)xSU(2) Wess-Zumino- 
Witten (WZW) theory. 

Our approach to analyse the scaling limit of the HLBA equations of the HF Hubbard 
chain is complementary to the existing ones |]IT?| , |TT| , |2"2"[ , and at the same time it also makes 
possible to reproduce most of the known results. For convenience we list here our main 
results: 

1. The resulting HLBA equations form three disjoint sets (for the massive excitations 
and the right moving ('r'), resp. the left moving (T) sector) each determining the 
parameters of only one kind of excitations. (This shows that only the same type of 
excitations interact in a nontrivial way.) All the limiting HLBA eqs. have the same 
structure: 

£ 21a,,- 1 ( ) = 2 T J x + £ 2tan-' ( , (1.4) 




7T 



i ln r U-*i) r ( 1 + *f 
1 r(| + zf)r(i-if) 



with p being the momentum of the excitations, x stands for the rapidity of the cor- 
responding kind of excitations and the As are parameters needed to distinguish the 
various internal symmetry states of the given kind of excitations. The actual quantum 
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numbers I x obey some "parity-prescriptions". For the massless 'r' and T particles 
these prescriptions depend on the number of particles in the given branch only, imply- 
ing that the massless sector decouples from the massive one in a sense. To be more 
precise the rapidities of the V (T) massless particles are independent of the massive 
sector and the state of the T (V) branch. (A prescription, however, which requires 
that the total number of physical particles must be even, still represents a nontrivial 
coupling between the sectors.) 

2. While the massive particles carry spin, the massless ones carry isospin (charge). We 
find, that for most states neither the V nor the T excitations form separately isospin 
eigenstates, only the 'r' and T sectors together. These states are characterized in 
addition to the parameters of the V and T sectors by an additional set of parameters, 
9, satisfying yet again BA-type equations: 

1 9 — 2 i $ + 2 x — , i 9 — 9 , . 

n r 2 tan -1 h n{2 tan" 1 = 2ttJ 9 + } 2 tan" 1 . (1.5) 

u u p 2u 

The 9 parameters do not appear in the equations determining the rapidities so they 
do not influence the energy, they are needed, however, to account for the isospin of 
the states. We argue that for these degenerate states, i.e. for those which differ in the 
number and values of the 9 variables only, a new basis can be constructed in which 
both the 'r' and '1' sectors are in isospin eigenstates. While for the original eigenstates 
the chiral charge is not well defined, the new basis consists of eigenstates of the chiral 
charge, i.e. the spectrum of the limiting theory is chiral invariant. 

3. By computing finite size corrections together with the fine structure of the spectrum 
in the scaling limit, the conformal properties of the limiting theory are exhibited. The 
massless sector is shown to correspond to a conformal field with central charge c = 1 
and Gaussian conformal weights 

A = -(n + m) 2 and A = — (ra — m) 2 , (1-6) 

where n and m are integers related to the numbers of the 'r' and '1' particles. Eq. ( |1.6|) 
lends further support to the identification of the massless sector with a level one WZW 
theory. We note that ( |1.6| ) together with the S-matrices for the massless sectors 
calculated in Section III. also yields convincing evidence that the massless bootstrap 
S'-matrix proposed for level one WZW models in Ref. [^TJ is correct. We also note, 
that the theory obtained in the scaling limit is an example of a conformal field coupled 
to a massive one. We find that the spectrum of the massless sector has a conformal 
tower structure even if the massive sector is not empty, but in this case the positions 
of the towers can be shifted. 

4. The behaviour of states with a finite density of excitations is also dicussed. The zero 
temperature free energies (which are nothing but the ground state energies in the 
presence of a chemical potential) of the different sectors are independent of each other, 
and the total free energy /(/i, v) can be written as 
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fO*,v) = -£*{ JL )--, (1-7) 

7T mo 7T 

where \i and z/ denote the chemical potentials for the massive resp. massless particles. 
We also exhibit the asymptotic series of $(/i/mo) (in l/[ln(/i/mo)]) for high particle 
densities. 

5. Some additional evidence is provided that the scaling limit of the HF Hubbard chain 
is the SU(2) CGN model (|Df) with g' = 0. First we show that when g' = the 
full symmetry of the CGN model is in fact SO (4) which is to be expected if it really 
corresponds to the scaling limit of the HF Hubbard chain. Secondly we show that 
the zero temperature free energy calculated in the high density limit by perturbation 



theory p3] is in complete agreement with the result obtained from the scaling limit 
of the HLBA equations. This also provides a rather strong (though indirect) evidence 
that the beta-function of the CGN model is correctly obtained (up to two loops) in 
the scaling limit. We have also made a one loop renormalization group analysis of the 
naive continuum limit of the HF Hubbard model which is neither a chiral nor a Lorentz 
invariant theory. Our analysis indicates that the renorm trajectory of the naive limit 
aproaches that of the CGN model, lending further support to their identification. The 
one loop result is, however, not conclusive as the renorm trajectories are driven into 
the strong coupling regime. 

The paper is organized as follows. First (Sec. [TT].) we review shortly the BA and HLBA 
equations describing the states in question. The relativistic limit of these equations together 
with their interpretation is given in the second part (Sec. pTI[), while the third part (Sec. |TV| .) 
is devoted to the discussion of the finite size corrections and conformal properties. In the 
fourth part (Sec. M.) some states with a finite density of excitations are discussed. In the last 
part (Sec. |VT|.) some arguments supporting the equivalence of the SU(2) (actually SO(4)) 
symmetric CGN model and the scaling limit of the half filled Hubbard chain are collected. 
At the end of the paper in Appendix |A] we illustrate the isospin structure of the massless 
states. In Appendix || we show that the approximations used to derive the HLBA equations 
for the discrete chain remain correct in the scaling limit too, and in Appendix |CJ the finite 
size corrections are discussed. In Appendix |D] we sketch the perturbative calculation of 
the free energy in the CGN model, while Appendix [E] contains the derivation of the naive 
continuum limit of HF Hubbard chain. Finally in Appendix [F] a one-loop renormalization 
group analysis of the theory obtained in the naive continuum limit is given. 



II. THE BETHE ANSATZ EQUATIONS 



Let us exhibit first the two commuting SU(2) symmetries of the Hubbard Hamiltonian ( |1 . 1|) 
explicitly P,|2~3||: the 'spin' SU(2) is generated by 



1 

2 
1 

2i 



c ll c j,d ' 



(2.1a) 



5 



J 

while the generators of the other SU(2) ('isospin') which is related to the charge are: 



i2 = iU-i) j (wu-4Ad, ( 2 - lb ) 

J 3= 9 EC 1 - «i,t - ■ 
Z J' 

We remark that the spin and the isospin can be mapped onto each other by the transfor- 
mation 

c j,l ~* %l ' S,T ~* ( _1 ) J . c Jt ' fo lcl 

c j,l ^ c j,l ' c j,T ^ ^ C i,t • 

(This transformation applied to the H changes the sign of U.) 

The diagonalizing the C P- -ID Hamiltonian is equivalent to solving the set of nonlinear 
equations 

Nk, = 2nL - Y 2 tan" 1 sin V~ A ° ( 2 . 2a ) 

V 2 tan' 1 Aa ~ Sinfcj = 2vr J a + V 2 tan" 1 . (2.2b) 

Here N e (< N) is the number of electrons and M(< N e /2) is the number of down spins, i.e. 
S z = (N e /2—M) and J 3 = (N—N e )/2. The solutions correspond to longest spin- and isospin- 
projection states (S 2 = S Z (S Z + 1), I 2 = h^h + 1)) provided all the A are (whatever large 
but) finite. The Ij and J a quantum numbers are integers or half-odd-integers depending on 
the parities of N e and M: Ij = M/2 (modi), J a = (N e + M+ l)/2 (modi). Once equations 
Q2.2| ) are solved the energy and the momentum of the corresponding state can be calculated: 

E = NU/A - Y^{2t cos kj + U/2) , P = ^ % , (2.3) 

3=1 3=1 



and also the wave-function can be given |25| 



For the considered U < attractive chain, near the ground-state most of the electrons 
form bound pairs with wavenumbers given (up to corrections exponentially small in the 
chain length (see Appendix [B|)) as 

sinfc ± = A±ra, (2.4) 

with u = \U\/4t as given earlier, and A being a subset of the set A. By this relation fc ± 
can be eliminated from Eqs. (|2.2a|) and (|2.2b|) and one finds that the k wavenumbers of the 
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unbound electrons, the As connected with the distribution of their (uncompensated) spins 
(i.e. those elements of the original A set which are not in the subset A), and the As of the 
bound pairs satisfy the equations M 



1-Klj = Nkj - 



E„ 1 Sill kn \fy - S -\ „ 1 Sill kn -A.77 /„ \ 

2 tan" 1 3 - - 2 tan" 1 3 - v - , (2.5a) 



n(\)+n(A) 
Ij = (mod 1) ; 



2 tan = 2tc J a + > 2 tan — - — - , 

u „-,2u 



j=i u 0=1 zu (2.5b) 

n (k) _ W (A) + 1 
J Q = (modi) ; 



27r J,, = (sin ^A^ - iu) + sin ^A,, + iu)j - 

n(fc) A — «in £• • n(A) A - A 

_ £ 2 tan" 1 " 5 J - J2 2 tan" 1 v u , (2.6) 
j=i u v=i 2m 

7l(fc)+7l(A) + l 

^ = ~ (mod 1) . 

Here n(k), n(X) and n(A) denote the number of unbound electrons, the number of unbound 
electrons with down spins and the number of bound pairs, respectively (N e = n(k) + 2n(A), 
M = n(X)+n(A)). (We display the prescriptions for the quantum numbers / and J together 
with the equations: in the following derivations the quantum numbers will be redefined by 
absorbing different constants into them and so the actual prescriptions will change from step 
to step.) The energy and momentum expressed in terms of these variables are 

E = -Nu - 2t (cos kj -u)- 



5^ 2t [Jl - (A r/ - iu) 2 + Jl- (Arj + iu) 2 - 2u 



V 

j ol r\ 



(2.7) 



These equations can be solved using methods by now standard in the treatment of BA 
type equations []?l |2"6l , [2"8f . The applicaton of these methods for the present case involve the 



following considerations and steps: 

i. In the ground state all of the N e = N (half filling) electrons are condensed into bound 
pairs, whose rapidities are given by Eq. fl2.6|) with n(k) = and J v being consecutive 
integers or half odd integers between ±(N/2 — l)/2. For a macroscopically large 
A^ the A,jS satisfying Eqs. (|2.6|) will be given by a distribution, cr (A), so that the 
number of the A^s within the interval (A, A + dA) is given by Nao(A)dA. Then the 
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summations over the As in (|2.6|) and in (|2.7|) for the energy can be replaced by integrals: 
X)„/(A„) — > j /(A)o"o(A)dA and the density, (Tq(A), is determined by a linear integral 
equation derived from (|2.6|). 



ii. There can be two typs of excitations. One is connected with breaking up pairs, that is 
introducing real ks and As. The other type is introduced by changing the distribution 
of As by leaving holes in the ground state distribution and introducing complex As. 
Also for the excited states the set of real As can be described by a density, cx(A), 
differing from the ground state density in 0(1/N) terms due to the excitations. After 
eliminating the real As from the set (|2.5| a-b) and ( |2.6| ) one is left with a set of equations 
for the unknown ks (which are real) and As, the positions of the holes in the real A 
distribution, A^, and the set of variables O n representing the complex As (Appendix 
[B]). The HLBA equations obtained this way read: 



- V ' A '.'~ / I + e -2H M J o(^) sin (^ sinA: j)~ J = 2?r/ J + 



-oo 



Ei sin kj — X a ^ , /sin fc,- — sin k~, 
2tan > J — -, 

« y \ 2u ) (2.8a) 

V2tan- 1 fth 7r(sin ^~ Aft) 

Ij = N-n{h)-( n (k)-2n(X)) ^ 



£ 2 tan- 1 Aa ~ Sinfej = 2vr J a + £ 2 tan" 1 , 

i , x " /s 2m (2.8b) 

n(fc) + n(A) + 1 
J Q = (mod 1) ; 



OC 



A / — smuA h — = 2irJ h + > 2 tan 

J 2chu;M uo „ 



u 



£ ( Ah ~ u Ah ' ) + £ 2 tan- (th * {Ah -™ k » ) , (2.9a) 



h> 



1 f N + n(k) + 2n(Q) - n(h) \, , 
J h — h 1 mod 

2 V 2 J v 



1) 



£ 2 tan' 1 9n ^ Ah = 2nJ n + £ 2 tan" 1 9ft Qn ' , 

* n(ft)+n(e) + l , '" (2 - 9b) 

J n = (mod 1) . 

Here n(/t) denotes the number of holes (i.e. the number of A/J and n(Q) is the number 
of variables 0. The energy and the momentum of the state given by the solution of Eqs. 
(pja-b) and (pfc-b) are 
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E = Ne + ]T e s {kj) + £e c (A h ) , 

3 h 



(2.10a) 



P = J2Ps(kj)+J2P^h) , (2.10b) 

3 h 



with 



£o = — 2r / tu , 

u; 

e a (k) = — It I (cos k — u) — I — ^ 6 cos(u; sin k)— \ , (2.11a) 

— oo ' 

, JAuj) . duj 

e c A = 2t / -^coscuA— , 

— oo 

and 




Ps(k) =\k— / : sm(a; sin k) — 

V J 2cna;M uj 

v — oo ' 



:2.11b) 



-oo 



(This result is the same as the one obtained from the repulsive case J7j using the comple- 
mentary solutions H of the Eqs. ( |2.2j a-b). The convention used for the indices s and c is 



that of Ref. [|K|.) The prescriptions for the I and J quantum numbers are expressed by the 
numbers of excitations only and numbers like n(A) are eliminated. This is possible through 
the relation 

n(h) = N - 2n(A) - n(Jfe) + 2n(6) , (2.12) 

which is a consequence of Eq. (|2.6|) . This relation also implies, that (for even N) the number 
of excitations n(k) + n(h) is always even. It also allows to express the spin and the isospin 
(charge) in a simple form 

S> = U{k) ~ 2n(A) , (S» = S'(S' + 1)), 



2 

n(h) - 2n(6) 



(2.13) 
{P = h{h + 1)) , 



which suggests, that the two kind of excitations carry either spin or isospin (charge) but 
not both. It is to be emphasized, that while both S z and Is are either integers or half-odd- 
integers, due to the above requirement their sum is always an integer. (After all this is a 
consequence of the fact that changing the number of electrons in the system by one will 
change both S z and J3 by plus or minus 1/2.) 
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III. THE RELATIVISTIC LIMIT 



The HLBA equations and the energy-momentum dispersions have the form of ( |2.8| a-b), 
( |2.9| a-b) and ( |2.11| a-b) in the half filled case only, i.e. using them automatically fixes the 
band filling. Since letting N — > oo means also N e —>■ oo, to keep the mass gap finite the 
interaction u has to — ► 0. In the u — > limit |JID|,PT| 



e s (k) 
Ps{k) 



(8t/ir)K 1 (n/2u)ch(n sin k/2u) 
(4/n)K (n/2u)sh(ir sin k/2u) 



(3.1a) 



and 



e c (A) = AtI x (>K/2u)e 



-tt\A\/2u 



Pc(A) 
(|A|>1) 



7T 



^nA{-- + 2I (n/2u)e-^/ 2u 



(3.1b) 



with Ko t i and Io,i being the modified Bessel functions. The lattice- or quasi-momenta p s 
and p c are dimensionless as they generate the discrete translations of a lattice. To obtain 
the corresponding momenta in the continuum, one has to divide them by the lattice spacing 
a. In the a — > limit, however, the ±7r/2a terms in p c /a would diverge. For this reason, 
before taking the a — > limit the lattice has to be redefined in such a way, that four lattice 
sites make one elementary cell. Then the quasi-momenta ±7r/2 are equivalent to zero, i.e. 
they can be dropped from the total momentum. Taking now the limit (|1 . 2j ) leads to 



E - Ne = £ e m {K) + £ e r {r,) + £ e,(0 

i 3 = EPm(«)+EPr(»7)+EPi(0 



(3.2a) 



with 



IA>1' 



Pr(v) 



lim 



lime s (/c) 
e r (?7) = lim e c (A)| 

ej(0 = lim £ C (A)| A< _ 15 p,(^) = lim 

e m(ri/) and p m(ri/) given by 

e m («) = m cosh(ft) , 
er(^) = ^exp(+?]) , 

e «(0 = -^ ex p(-£)> 



limp s (fc)/a, 
p c (A) + vr/2 



Pc(A) - vr/2 



A>1 



A<-1 



7rsin k 

2u ' 
tt(A-2) 

2m 

tt(A + 2)^ 
2u : 



(3.2b) 



p m («) = m sinh(/t) , 

MO = — ^ ex p(-0 • 



(3.3) 



limit the energies and momenta are 
-> +or-2. The limits of the HLBA 



As one can see from ( |3.3| ) and ( |3.2| a-b) in the u — > 
finite, if k, rj and ^ stay finite, i.e. k —> and A 
equations are compatible with these requirements. 

The HLBA equations for the massive particles are obtained from Eqs. (|2.8a ) and ( 2.8b ) 
by substituting (|2.11b| ), ( p.2b|) and \ a = n\ a j2u into them: 
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a 7T/2 y V 7T 



The last two terms give a constant 



lim ^2 tan" 1 (th 




-(n(0 -71(77)) 



(3.5) 



with 77(77) resp. n(£) denoting the number of 77s resp. £s. The constant ( |3.5| ) can be absorbed 
into Jj (from which N/4 can be dropped, as due to the redefinition of the lattice, N must 
be a multiple of 4). Finally using that 72(77) + n {0 = n(h) and renaming n(k) and 77(A) to 
ti(k) and n(x) respectively, one obtains 



L Pm (^) = 2,1, + £ 2 tan" 1 - £ {^ £ ) 



J. = ± W (^)+ W (0 + K«)- 2W (X)) nfo) - 77(0 (mod x) (3 . 6a) 

J 4 4 

n r + 77; + (n(«) - 2n(x)) . n r -ni 

H (mod 1) ; 



4 4 
(with n ri / defined later in ( |3.9a| ) and ( |3.10a|) ), and 

y 2 tan" 1 = 2tt J a + y 2 tan" 1 Xa - X(3 

7 ^ 7 

J a = " (/t)+n(x) + 1 (modl). 



(3.6b) 



We have written Ij in this form (using that the first term is either an integer or a half- 
odd-integer) to emphasize its symmetry with respect to 77(77) and n(£) (n r and 77;). This is 
actually the left-right symmetry of the system: interchanging 77(77) and and changing 
the sign of all k and x yields an admissible solution (i.e. with proper quantum numbers). 

The limit of the HLBA equations for the massless particles is constructed in a similar 
way. First we observe that as u — > Eqs. ( |2.9a| ) force the A/jS to take their values around 2 
and —2. At the same time some of the condense around 2 and —2, but some may have 
values |G| qk 2. So the set of 6 can be replaced by three sets: 

= -tt(0 - 2)/2u (6~2), 

cp = -n\e + 2)/2u (6~-2), (3.7) 
9 = Q (lim(2- |0|) ^ 0) . 



(Their numbers satisfy 77(7?) +n(cp) + n{6) = 77(G)). In the variables (|3.7|), equations ( |2.9a| ) 



and (|2.9b| ) split into three disjunct sets which can be treated in the same way as the set 
describing the massive particles. For the right-going particles (77 and 7?) one has 
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N(l - «*(,,)) = 2,J f - E2t»-' 2^ + X>(^) 
£ 2 tan- f 5£J1^ _ £ 2 tan -i I H - 2 ~ «* 

k 



7r/2 u) „ \7r/2 U 

+ EH^-D-E2Un-( t h( 



4m 2m/ 



(3.8a) 



T 9 1 ^ + »(«) + - n(g) + 2n(y) - n(fl + 2n(6>) \ 
f = 2 [ 2 j t modl J ; 



_ V 2 tan" 1 - V 2 tan" 1 f - -) = 2nJ q - 

y 7T/2 y \ 71 uj 9 

V2tan- 1 ^-^-V2tan- 1 ^- V?fc -^ 



- e 2 tan ~ 



7T 2"U 



J 3 = (mod 1) . 

By absorbing all those terms into J/ and J g , which become constants in the u — > limit 
(and using that ti(k) + n(£) + 77(77) = wen) leads to 



L Pr (rj f ) = 2nJ f + £ 2 tan' 1 - £ 

7T/2 „ V 7T 



» /' v 7 (3.9a) 

/77 r 1\ 

J/ = - f -j + - J (mod 1) , n r = 77(7/) - 2n(7?) ; 
£ 2 tan" 1 ^-Z^ = 2tt J g + E 2 tan' 1 tl^llL , 

/ 7r/2 / ^' n (3.9b) 

7777+7779+I 
J g = (mod 1) . 

An analogous procedure yields for the left-going particles (£ and ip) 
L Pl (Q = 2nJ q + £ 2 tan" 1 - E > 

E 2 tan- 1 ^-^ = 2. J k + E 2 tan" 1 , 

9 7T/2 fc , 7T 

J fe = (mod 1) . 

Finally the variables 9 are defined by the u — > limit of the solution for the system 



(3.10a) 



(3.10b) 
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n(r/)2tan" 1 — - n(?9)2tan" 1 — -+ 

u 2u 

n U) 2 tan' 1 - n ((p)2 tan" 1 = 2rt J n + V 2 tan' 1 ° n - ° n ' , (3.11) 

it 2m , 2m 

7^(77)+ n(£)+n(tf)+n(</?)+n(0) + l 
J„ = (mod 1) . 



As for 9 96 ±2 



2 tan" 1 = 4 tan -1 =f tt + 0(m 2 ) (mod 2tt) , (3.12) 

2u M 



the m — > limit of solutions of ( p.ll|) is the same as the u — > limit of solutions for the 
equation 

M r 2 tan -1 — h M/2 tan -1 n — = 2i\ J n + 2 tan -1 — — , 

u n , 2u (313) 

n,. + n , + n (0) + 1 
J n = (mod 1) . 



We remark that Eqs. (|3.13|) have the same structure as the BA equation of the chiral invariant 



Gross- Neveu model [18]. We also note that the 9s do not influence the energy (not even 
indirectly through some other quantity). They cannot be discarded, however, as they are 
important to give the isospin-state corresponding to a solution. 



Equations (|3.6|a-b), (|379|a-b) , ( |3.10| a-b) and ( |3.13| ) are the scaling limits of Eqs. (|2T8|a-b) 



and (|2.9| a-b), and they are the secular equations of the limiting model. Their solutions are 
meaningful under certain conditions. Firstly all roots have to be different: if two ks, As 
or two $s, 9s or ips are equal, the wavefunction is identically zero, thus such solutions are 
(although formally possible) meaningless. Similarly no two r/s or £s can be equal as these 
variables are originally holes in the A distribution. Based on the fact that (|3.6b| ), ( |3.9b| ), 



( |3.10b|) and (|3.13|) are very similar in structure to the original equations ( |2.2b| ), we expect 



them to have meaningful solution if n(r]) — 2n($) = n r > 0, n(^) — 2n((p) = m > and 
n r + rii — 2n(6) > 0. In addition to these conditions 2 — \8\ ^ should also be required for 
all 9 (as this has been used in deriving ( |3.13| ) (any = ±2 should be represented by a ^ or 
if)). Constructing the simplest solutions of ( |3.13|) we have found that this requirement can 
be satisfied if n(6) < min(n r ,nj). 

We interpret Eqs. ( |3.6| a-b), ( |3.9| a-b) and (|3.10| a-b) as the secular equations of the massive 
('m') and the two (V and T) kinds of massless particles, and we argue that Eq. (|3.13|) 



describes how the two massless sectors combine. In connection with their structure and 
content we note the following: 

i. Eqs. ( |3.6| a-b) are of the BA form: this structure is typical for the HLBA equations 
of systems with excitations having an internal degree of fredom. In the case of the 
massive particles this degree of freedom is the spin (see ( |2.13p ). A solution of the 
secular equations corresponds to a state 

S 2 = l(l + 1), S z = m, I = m = n{K) ~ 2n{x \ (3.14) 
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i.e. there are ti(k) particles, n(n) — n(x) with up and n(x) with down spins. States 
with m < I can be constructed from these by acting on them by the operator S~ oc 
(S x — iS y ). If this operator is applied v times, we still have I = (n(n) — 2n(x))/2 but 
m = I — v. This way, if ( |3.6| a-b) give all the / = m states, all the \m\ <l states can be 
constructed. These states are degenerate with the corresponding I = m one. It should 
be noted, that the action of the S~ operator can be implemented in the BA formalism 
by introducing as many extra x s as many times the S~ acts and taking these extra x s 
to + or — oo. 



j|a-b), ( |3.10| a-b) and ( |3.13|) describe the massless sectors where the internal degree 



of freedom is the isospin, with its third component being the charge (particles with 
J 3 = ±1/2 have charges ±1). Since the secular equations of the two massless (T and 
'r') sectors are independent of each other (see item vi. too) and have the same structure 
as those of the massive particles, it is natural to think that the isospins carried by the 
'1' and 'r' particles (based on the analogy to the spin of the massive sector) can be 
defined separately. On the other hand we can make exact statements only about the 
total isospin (( p,16b| ), ( |3.17b[ ), ( ^.18b| ) and ( |3.19b| )), but not about the isospins of the 
'1' and 'r' sectors separately. A reason for this is that in the original Hubbard model 
the '1' and V sectors (so their isospins) are not even defined, as the T and 'r' sectors 
separate in the scaling limit only. In the following we try to resolve this apparent 
contradiction and propose an interpretation according to which it is possible to choose 
a basis, where both the T and 'r' sectors are in isospin eigenstates decoupled from 
each other. 

A central element of our proposed interpretation is that in any state characterized by 
n r and ni both massless sectors are in eigenstates of the (J) 2 with eigenvalues 

( /W ) 2 = T (t + and = + (3 ' 15) 

This is supported by the following. If one of the sectors is empty (the V one, say), the 
other (the T sector) carries the total isospin and its state is indeed an eigenstate of 
(J^) 2 and I3 , where the value of (J®) 2 is the one given above and I® = rii/2. Acting 
on this state by I oc I\ — il<i decreases 7g by one, leaving (J^) 2 (and all the other 
quantities) unchanged^ By a repeated application of I~, all the I-Z3 | < rii/2 states 
can be constructed, i.e. a single solution of Eq. ( |3.10| a-b) represents a complete set 
of degenerate states with (7^) 2 given above and I-Z3 | < rii/2. Analogous statements 
hold for the V sector if the '1' one is empty. Based on the fact, that the solutions of 
( |3.9| a-b) and (|3.10| a-b) are independent of each other, it is very plausible to assume that 



this also holds for the states corresponding to them. Thus we propose that any state 
corresponding to a solution of Eqs. (|3.9|a-b), ( |3.10| a-b) is built up as a combination 



As in the scaling limit / may be not well defined, to apply I means in fact to act on a state 
of the finite chain by I~ first and take the scaling limit afterwards. 
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of the 'r' and T states represented (in the above sense) by the given solutions of 
Eqs. (|3.9|a-b) resp. (|3.10| a-b), even if none of the two sectors is empty. 

The states with n(9) = satisfy, in addition to Eqs. ( 3.15 ). 

= 4" = |, (3.16a) 



consistently with the exact formula for the total isospin, 



Acting on such a state by I~ decreases /3 by one, leaving the value of I 2 unchanged. 
The resulting state is a symmetric linear combination of two states in which both 
(I^) 2 and (J^) 2 are unchanged, but in one I3 , in the other I3 is decreased by 1. 
One can act with I~ on a state several (n r + n{) times. If I~ is applied v times, the 
resulting state is given as such a linear combination of the states with 

r (r) n r n r \ m Til »A , v 

J 3 = y - ^ ^> -y J , 7 3 - y - - yj > "r + *i = " , (3.17a) 

for which the total isospin is 

2 n r + nifn r +ni \ n r + ni - 2u ,„.„,, 
— + 1J> ^3 = g • (3.17b) 

If n{9) 7^ 0, none of the sectors are in eigenstates, but such a state is built up as 
a suitable linear combination of states characterized by ( |3.15j ) and 

4 r) = y-^, 4° = |-^, i/ r + i4 = n(0), (3.18a) 

so that the total isospin is 

2 n r + ni - 2n(9) f n r + rii - 2n(9) \ n r + m - 2n(0) /0101 x 
/= g ^ g + 1J, /a = 3 , (3.18b) 

as required by ( p. 13 ). (This interpretation is consistent with the following: a 9 can 
be considered as a 1? ~ 00 or as a 1/3 ~ -00, and (in analogy to what is said in item 
i.) the presence of such infinite variables can be identified with the action of 
resp. .) Of course, if n r + n\ — 2n{9) > 0, I~ can also be applied to these states 

further decreasing J 3 while leaving J 2 unchanged. This way, the action of I~ v times 
(n r + ri[ — 2n{9) — 2v > 0) results in a state, which is such a linear combination of the 
states with 

4 r) = y-^r, 4° = y-^, i/ r + 14 = n(0) + 1/, (3.19a) 

for which the total isospin is 
p _ n r + m- 2n(9) f n r + m - 2n(9) + \ ^ _ n r + n { - 2n(9) - 2v ^ 
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All the states with different n(6) and v corresponding to a given pair of solutions of 
Eqs. (|3.9|a-b) and (|3.10| a-b) are degenerate, but linearly independent. Based on the 



above interpretation, it is possible to choose such a basis in the space of these states, 
in which the basis vectors are labelled by the quantities (I^) 2 , I3 , (i = l,r) (see also 
Appendix |A|). We also note that in this basis the length of the total isospin 
is not a good quantum number any more. 

iii. The chiral charge is given by C = 2(1$ — I®). For a large class of BA eigenstates 
(described by fl3.17| ), ( |3.18| ) and ( J3.19Q ) the chiral charge is not well defined, as these 
states are not eigenstates of C. This is a consequence of the fact, that in the lattice 
model the chiral charge is not even defined. The new basis vectors described in item 
ii. are, however, not only eigenstates of the energy and momentum, but also of the 
isospins of the two sectors separately. For this reason they are eigenstates of the chiral 
charge too, ie. the spectrum is chiral invariant. The possible eigenvalues of the chiral 
charge C for states characterised by n r and n\ are 

C = n r + ni — 2n , n — 0, 1, 2, • • • , n r + rii . (3.20) 



iv. The scattering phaseshifts of the physical particles can be reconstructed from their 
secular equations. The method for this is based on the idea, that the deviations of 
the particle momenta from the free values can be interpreted as the phaseshifts of 
the particles scattering on each other [EfflE7lEBl. Consider a two particle scattering- 



state on a ring. If the momenta are p\ and P2 then periodic boundary condition 
requires Lp\ + 5\2 = 2nni and Lp 2 — 5\2 = 2im 2 with n% and n 2 integers and 5 12 being 
the phaseshift. Writing the scaling limits of the HLBA equations in this form the 
phaseshifts can be found. 

A triplet state of two massive particles is described by two ks and no x s - For such a 
state Eq. fl3.6a| ) yields 

5 12 = <f> (f) + , X = K 1 -K 2 , (3.21) 

where the 7r comes from the parity prescription for the parameters Ij. A singlet state is 
characterised in addition to K\ and K2 by a x f° r which Eq. ( |3.6b| ) yields x — (/ci+/C2)/2. 
For this case the above reasoning leads to 



'12 



' ' 2 tan" 1 (-) . (3.22) 



7T/ \7T 



^From the phaseshifts ( 3.2 1|) and ( 3.22|) the S'-matrix of an 'm-m' scattering can be 
given as 

- , N f. , fx \ 1 xl — iirll ln rio . 

S mm {x) = - exp \i(p - )) ; — , (3.23) 

I \7T/ J X — 171 

where I and ii are the identity resp. permutation operators acting on the spins of the 
two particles: 
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r'i a 2 — S A T\< a 2 = A A (3 04) 

In a similar way the S'-matrices for the 'r-r' and '1-1' scatterings are found to be 

S„(x) = S u (x) = -S mm (x) (3.25) 

with I and II acting on the isospins, and x = rji — r]2 resp. x = £i — £2 for the 'r' resp. 
T sector. The scattering matrices ( |3.23| ), ( |3.25| ) coincide - up to an overall sign - with 



those given in Ref. ||10|| . (Actually it is argued in Ref. that in general S'-matrices 
can be determined only up to an overall phase from the HLBA equations, and further 
considerations are needed to determine this phase ambiguity.) 

v. The 'm-r', 'm-1' and 'r-1' phaseshifts do not appear in the equations explicitly. Since 
they are just constants, we absorbe them into the quantum numbers Ij, Jj and J q . 
The prescriptions for these quantum numbers are compatible with these phaseshifts 



being + or — 7r/2 as obtained in Ref. [10 



vi. We note, that the 'r' and '1' massless sectors are apparently independent of each 
other and the massive sector in the sense that in the Eqs. (|3.9|a-b) and ( |3.10| a-b) 
determining the rapidities of these sectors no parameters of the other sectors appear, 
i.e. the rapidities rj and £ are independent of the states of the other sectors. This 
is a consequence of the requirement that the total number of excitations should be 
even (Eq. (|2.12| )). To see this consider a left moving (T) particle. In the equation 
determining its rapidity (besides the phaseshifts due to scattering on the other '1' 
particles) the sum of the phaseshifts corresponding to scattering on the massive and 
£ r' particles appear. This sum can be expressed (mod 2n) by the number of the '1' 
particles only, and this is why no parameters referring to the 'm' and 'r' sectors appear 
in ( |3.10| a-b). Similar statements hold for the 'r' sector, but not for the equations of the 
massive particles. The phaseshifts of the 'm-1' and 'm-r' scattering have an opposite 
sign, so it is the difference of the number of the '1' and 'r' particles which enters into 
Ij of ( |3.6a|) . This suggests that in spite of the apparent independence of the massless 



sectors the different sectors are nontrivially coupled. (This is also indicated by the 
prescription that the total number of excitations should be even.) 



IV. CONFORMAL PROPERTIES 

The existence of a massless sector signals an underlying conformal field theory. By 
calculating the finite size corrections to the ground-state energy and the fine structure of 
the spectrum arising due to the quantization of the momenta (tower structure) one finds 
the central charge and the conformal weights [[H]]. To find the finite size corrections to 



the ground-state (dressed vacuum) energy to the spectrum one has to calculate them for 
the Hubbard chain and take the scaling limit afterwards. In this Section we give first the 
finite size corrections to the spectrum of Eq. (|2.10a|) , then we take the limit ( p..2|) and finally 



we show that the resulting spectrum has the conformal tower structure. (An alternative 
treatment is to consider the lowest energy states in the finite N Hubbard chain, and take 
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the scaling limit of these states directly (as it is described in Appendix Both treatments 
we present are based on an Euler-Maclaurin type summation formula, but there exist other 



methods too I02LI33 



When deriving equations (|2.8| a-b) and ([2.9|a-b) and the spectrum (|2.10| a-b) from 
Eqs. ( |2.5| a-b) and ( p.6| ), the summations over the As have been replaced by integrals: 

(l/N)J2 v f{\) J/( A ) (7 ( A ) rfA witl1 bein g the density of the A. This replacement 
becomes exact in the N — > oo limit only. For large but finite N using an Euler-Maclaurin 
type formula even the (1/A) fi-h-r)) — I /(A)o"(A)dA difference can be taken into account. 
The corrections obtained this way are referred to as finite size corrections. The method is 
elaborated to quite an extent JftJ, and it can be adapted for the present treatment (details 
will appear elsewhere). As a result (|2.10a|) gets modified, an additional term appears: 



and e c (Ah) differs from e c (A^) by terms 0(1/N). In the limit (|1.2| ) t can be replaced by 
l/2a, and after taking u —>■ we have 

E-Ne = + J2 + E £r(v) + E 5(0 , 

QL (4.2) 

p = EM*) + Eft-fa) +Ea(fl • 

For energies of the order of unity e n / = e r> i (and p r j = p r ,i), but for small energies (of 
the order of 1/L) the form of these quantities slightly deviates from those valid for larger 
energies. It is important to note, that 

e r = p r and l\ = -pi (4.3) 

still hold. The secular equations for the massive particles is not effected by the finite size 
corrections, but those for the massles particles are. For the simplest case (when there are 
no i?s and (ps, these states are of the lowest energies within a tower) they read 



LPrfrf) = 2nJ f ~ E (Jlf, Vf>) ■ 

f 

J / = -(^ + I) (modi), 

and 

Lpi{Q = 2nJ q - E ^ (£s> &) > 



(4.4) 



(4.5) 



Also here, for energies of the order of unity <ft(x,x') = <p(x,x'), but for small energies the 
two quantities deviate (but even for small energies (x, x') = —(f> (x', x)). We note, that the 
structure of these equations, and that of the energy and momentum are the same as they 
would be without the finite size corrections, and to calculate the spectrum of the simple 
states we do not even need the actual form of the e r i, p r i and (j). 
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If there are no excitations, on the right hand side of the first equation in ( |4.2|) only the 
first term remains, showing that the universal part of the ground state energy of the massless 
sector is the same as that of a Gaussian field with conformal anomaly c = 1. 

Next we argue that the contribution of the massless excitations in (|4.2|) has the conformal 
tower structure, provided the variables r\ and £ satisfy Eqs. ( |4.4| ) and (p~5|). By summing 
these equations we get 

£^) + £*(0 = ^W. + Ji), 

with 



^ = £<//, Jl = ~Y.h- (4-7) 

It is not hard to see (based on the prescriptions for the J quantum numbers), that 

n r (n r + 2) ni{m + 2) 
J r = (mod 1) , Ji = (mod 1) . (4.8) 

Consequently the states with given n r and ri[ form towers: at a fixed n r all possible J r (n r ) 
differ by integers, i.e. J r {n T ) = J r (min)(n r ) + N r and similar statements hold for Ji(rii). Thus 



2tt 2tt 

-j^(Jr + Jl) = -j-(Jr(min)(n r ) + Jl( m i n )(ni) + N r + N t ) 

2tt 2-k y 

~j^{Jr ~ Jl) — -j-{Jr{min){n r ) ~ Jl{min){ni) + N r — N[) 



(4.9) 



J r takes its smallest value, if 



T TXf 1 fly 3 Tlr- 2Tir 1 / , _ „ \ 

J > ] "T + 2- -T + 2'"'--T + — • (4 - 10) 



as if Jf( m in) < —n r /A + 1/2 Eq. ( |4.4j) has no solution. Then 

Jr(min)(n r ) = . (4.11a) 

In a similar way we have 

2 

Jl(min)(ni) = y , (4.11b) 

thus for a given n r and rii the minimal energy AE and the corresponding momentum AP 
of the massless excitations (the apex of the {n r ,rii} tower) are 

According to these the conformal weights are 
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i.e. the corresponding scaling dimensions and spins are given as 



nt 



x 



+ 



nt 



nz 



nt 



4 4' 4 4 

(Here we use the notations of Cardy [31].) Introducing 



n r + rii 
n = , m 



(4.14) 



(4.15) 



Eqs. ( |4.14j ) takes the form 



x 



- (n 2 + m 2 ) , 



nm . 



(4.16) 



This structure formally coincides with the expected one for a c = 1 conformal field theory 
coresponding to an SU(2) WZW theory of level one, however, this correspondence is not 
complete since n and m can also take half-integer values! When the number of massive 
particles is even both n and m are integers, while if the number of massive particles is odd, 
n r ±ni is odd, and both n and m are half-integers. (This is a consequence of the requirement, 
that the total number of excitations must be even, as given by ( |2.12j ) and is also discussed in 
item vi. of the previous section.) Thus we can interpret ( |4.16|) as follows: when the massive 
sector is empty the massless sector can be described by an SU(2) WZW theory of level one. 
When the massive sector is not empty one cannot expect a simple conformal field theoretic 
description because of the coupling of the massless particles to the massive ones (as discussed 
in items v. and vi. of the previous section). Also in this case the energy contribution of the 
massless sector shows a tower structure, but the interpretation of the apexes of the towers 
as scaling dimensions and spins is not straightforward. 

It is to be noted, that in the states with no 8s the parameters n and m can be identified 
as n = I3 and m = C/2 with C being the chiral charge. For states with a finite number of 
9s such an interpretation is not possible. 



V. STATES OF FINITE DENSITIES 

Although in our derivation it has been supposed, that the number of excitations is small, 
the resulting equations (|3.6| ), (|3.9| ) and (|3.10f) can be used to describe states with a finite 
density of excitations too: even if the number of excitations is macroscopic, i.e. ~ pL, that 
is small compared to N since pL/N = pa — > in the scaling limit. In this section we discuss 
the lowest energy finite density states. As the different sectors are 'almost decoupled' (Sec. 
p|.), this can be done for the different sectors separately. (The coupling made explicit by 
the parity-prescription for Ij can give a contribution 0(1/L) in the energy density.) 

It is plausible to assume that for the lowest energy state of the massive sector (at a fixed 
number of particles) the Ij quantum numbers are consecutive integers or half-odd-integers 
centered around the origin. If in such a state the number of ks is macroscopic (N m ), they 
can be described by a density <j(k) determined from Eq. fl3.6a|) as 
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It is then clear, that the density will be the highest possible around the small ks (small 
energies) if there are no x s - (Strictly speaking this argument holds for real x s on ly ; but it 
can be modified to hold for more complicated \ configurations too.) After replacing the 
sum over the Kj>s by an integral we find that (t(k) satisfys the equation 

B 

(t(k) — / drf 'K(k — k 1 ) ct(k') — — -cosh/?, (5.2) 
J 2n 

-B 

where the kernel K(k) can be written in terms of the <fi of ( |1.4j) as 
The ct(k) is also subject to the condition 

N m 

dKa(n) = — — = p. (5.4) 

Lj 

-B 

In terms of the solution of Eq. ( |5.2|) the energy density £(p) is given as 

B 

£{p) = m J dn <j(k) cosh k . (5.5) 

—B 

(Equations ( |5.2[ - |5.4|) are exactly the same as the ones obtained in Ref. ]29| directly from the 
BA equations of the CGN model.) We present the explicit solution of Eqs. ( |5.2| - |5^ ) for a 
large density (p/m ^> 1) using the results of Refs. f23|j34| . To apply those results directly 
we also give the integral equation determining the zero temperature free energy, obtained 
from Eqs. (5T2) and (|5.4|) by a Legendre transformation: 

f(p) :=min (£(p) - pp) , (5.6) 
where p denotes the chemical potential. After a straightforward manipulation one finds that 

B 

f(p) = / dn e(/?) cosh k , (5.7) 



-B 



2tt 

where the function e(n) and the 'Fermi point', B, are determined by the equation 



[k) — J dn'K^K — k) e(/t') = p — m cosh k , (5.8) 



21 



together with the condition e(±B) = 0. It is now a simple matter to apply the results of 
Ref. [^,0] to find the solution for a large density of Eqs. ( |5.8|) . The asymptotic expansion 



of free energy is found to be given as: 



m = -f • 2 

ZTl 



1 



In ln(p/m ) 



.4 



21n(p/m ) 4 In (p/m ) 8 In (p/m ) 



In ln(/i/mo) 
ln 3 (p/m ) 



(5.9) 



where the constant A is 



A= l + 21nvr-61n2. 



(5.10) 



The energy density of the system is easily derived from (|5.9| ) using the definition of the free 
energy with the result 



Tip 



1 



lnln(p/m ) + 3 - 21n(2vr) 



21n(p/m ) 4\n 2 (p/m ) 8 In" (p/m ) 



lnln(p/mo) 
In 3 (p/m ) _ 



(5.11) 



The leading term agrees with that of Ref. the non leading terms disagree, however. 

The treatment of the massless sector due to the linear dispersion is much simpler than 
that of the massive one as we need not solve an integral equation. Consider for defmiteness 
the 'r' sector (the treatment of the T sector is completely analogous). In the lowest energy 
states there are no ^s, so if the number of particles is N r = g r L the J/ quantum numbers 
are (just as in ( [4.10| )) 



Jf 



N r 1 

T + 2- 



N r 3 



N r 2N r 



(5.12) 



leading to the energy density of the V sector S r (g r 
energy density of the massless sector is 



ttqI/2 and if g r = Qi = £?/2, the total 



S(g) 



Tig 



(5.13) 



The free energy density of the massless sector (defined as the Legendre transformation of 
S(g) with a chemical potential v) is 



fly) = - 



71 



(5.14) 



We note here that the above results can also be obtained directly from the Hubbard 
chain, by constructing a state with a finite density of excitations in the presence of chemical 
potential and magnetic field and take the relativistic limit afterwards. 



VI. RELATION TO THE CHIRAL GROSS NEVEU MODEL 

In this section we present some evidence supporting the equivalence of the theory defined 
by the scaling limit of the Hubbard model with the SU(2) CGN model. In fact, based on 
the perturbative and non perturbative results it seems to as that the the full renormalized 
CGN theory can be defined as the scaling limit of the HF Hubbard chain. 
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The Lagrangian of the SU(n) symmetric CGN model 
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can be written as: 



£ = i4> a 04>a + ~g S (4>a4>a) 2 + {^alb^af 



(6.1) 



where the ip a are Dirac fermions, a = 1 , . . . n and summation over the repeated a in- 
dices is understood. The theory given by the Lagrangian fl6.1| ) is invariant under a global 
U(n)®U(l) c hirai transformation; 



4>a 



(6.2) 



with M GU(ra) and a being a (real) constant. (We note that the couplings (g, g') in flO] 
are related to those of (|6.1| ) as g = g s , g' = g v + g s /2.) 

As mentioned already in the Introduction the CGN model is asymptotically free ||13 



and as first found in the large n expansion its spectrum contains massive particles |T2| 



Furthermore due to an infinite number of conservation laws preventing particle production, 
the S-matrix of the CGN model is factorizable, i.e. the multiparticle S-matrix can be 
written as a product of the two-particle S-matrices. Using the 'bootstrap' approach (when 
the spectrum of the theory is postulated) the S'-matrix of the model had been proposed 
long time ago [14|-|16|. From now on we concentrate on the SU(2) symmetric CGN model, 
whose massive spectrum is particularly simple as it consists of a single doublet. Note that 
in this case the particles are their own antiparticles. 

The bootstrap approach led to the following two-particle S'-matrix of the massive doublet 
in the symmetric channel Jy| : 



S(x) 



T(l + zx/27r)r(l/2 - ix/2n) 
T(l - ix/2vr)r(l/2 + ix/2n) 



(6.3) 



where x denotes the rapidity of the particles. (Note that S(x) coincides with exp{i<f>(x / V)} 
where <p(x) is the 'phaseshift'-function in Eq. ( |1.4| ).) 

As it has been argued in Ref. [0] using bosonization there is no contradiction between 
the existence of massive particles and the necessarily unbroken (continuous) chiral symmetry 
of the theory. The following picture of the CGN model has by now emerged [ 18 , 36 1 : the 



theory decouples into a massless and a massive sector. The massless excitations carry a 
U(l) charge but are singlets under SU(2). The massive sector consists of particles with 
a dynamically generated mass transforming non-trivially under SU(2) but without a U(l) 
charge. 



The application of the BA to the Hamiltonian of (|6\T 



18|JT9|| was a very important 
development as it made possible to obtain exact, nonperturbative results. In Ref. [0 the 
theory has been regularized by filling the Dirac sea up to a certain depth K (in a finite 
volume). This regularization is rather unorthodox as compared to other schemes, and it 
would be clearly desirable to work with a more physical cutoff. The main problem with 
other cutoff schemes is that they spoil integrability in general. Therefore it appears to be of 
particular interest if one could find an integrable lattice regularization of the CGN model. At 
first sight it is not immediately obvious that the Hubbard model would be a good candidate 
for a lattice version of the CGN model. The naive continuum limit (a — > 0) of the half filled 
Hubbard model is given by the following Lagrangian (See Appendix [E]) 
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(6.4) 



which significantly differs from ( |6.1| ), as the theory corresponding to ( |6.4| ) lacks chiral in- 
variance and it is not even Lorentz invariant. The theory defined by (|6.4|) is in fact rather 
similar to the (nonchiral) Gross-Neveu model. The lack of chiral invariance in the naive 
limit is not surprising as in the half filled case the umklapp processes play an important 
role. Nevertheless there is a large body of evidence that the scaling limit of the half filled 
Hubbard model is actually the SU(2) CGN model. We list some of them here and then come 
back to the problem of the naive continuum limit. 



The present analysis of the HLBA reproduces the expected spectrum of the CGN model 
consisting of a massive doublet and two (V and '1') massless doublets. The calculated 
phaseshifts are in accordance with the previously obtained results for the massive sector 
|L5| ) pO[| and with those of Ref. [|K|] where the complete S matrix (including the massless 



! r' and '1' doublets) has been found. The results show that the massive and massless 
sectors are practically (but not completely) decoupled. The massless sector is far from 
being trivial as the left- resp. right-movers scatter on each other exactly in the same 
way as the massless scattering theory proposed for a level one WZW theory |21| . The S 



matrix shows an SU(2) xSU(2) (in fact SO(4)) symmetry which is in fact to be expected 
as this is a symmetry of the Hubbard model which is preserved in the scaling limit. 
It is less clear, how such an SO (4) symmetric S matrix and spectrum can correspond 
to an SU(2) symmetric model. In fact this discrepancy is removed by the observation 
that the SU(2) CGN model also possesses a (hidden) SO (4) symmetry. To exhibit this 
hidden SO (4) symmetry we introduce four Majorana fermions x% (i = 1 • • • 4) instead 
of the complex SU(2) doublet ip a (a — 1 , 2) ipi = Xi + ^2 — X3 + Then using 
the Majorana properties together with the Fierz identities it is easy to show that there 
are two inequivalent SO (4) symmetric four fermion interactions in two dimensions: 

A (0) = fcx,) 2 4 1} = e^xasXjXkKXl • (6.5) 



In terms of Dirac fermions these interactions can be written as: 



(0) 



{Ma) 



£ 



(1) 



(lpal5lpa) 2 + -(^al^a) 2 



(6.6) 



Comparing ( |6.1|) and ( |6.6|) one can see, that the SU(2) CGN model with g v + g s /2 = 
possesses a hidden SO (4) symmetry (in analogy to the nonchiral SU(2) Gross-Neveu 
model whose full symmetry is in fact 0(4)). It is worth mentioning at this point that 
on the other hand in the scaling limit of HF Hubbard chain, eigenstates of the chiral 
charge can be defined (i.e. in that limit the chiral symmetry not present in the lattice 
model, appears). 

• The detailed analysis of the finite size corrections in section IV shows in a rather 
convincing way that the massless sector corresponds to a c = 1 CFT precisely at the 
SU(2)xSU(2) symmetry point, corresponding to a level one WZW theory. 
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The calculation of the free energy for macroscopic densities shows again the decoupling 
of the massive and massless sectors. More importantly as the free energy is a physical 
quantity it can be directly compared with the results of a perturbative calculation in the 
CGN model p3f valid for high densities because of asymptotic freedom of the model. 
Using the results of Ref. |23[ the free energy calculated in perturbation theory can 
be easily shown to be in complete agreement with that of the HLBA equations Q5.9D , 
( [5.141 ). In Appendix [B| we outline the strategy of the calculation of the perturbative 
free energy and simply state the result. 

Comparing the free energy ( |5.9j ) obtained from the HLBA equations with the pertur- 
bative one ( pi 0| ) one has to express the chemical potential /i in terms of h and Q in 
Eq. (|5.9|) . Since /i = hQ/2 it is easy to verify that the perturbative result (|D10|) is in 
perfect agreement with the exact one (|5.9|) provided 



m 



1 



(6.7) 



• The beta-function of the coupling u of the theory defined by the scaling limit of the 
HLBA system of the half filled Hubbard model can be immediately read off from 
Eq. (O): 

/?(«) : = ^ u = -/#V - pfV + ... fx=-, (6.8) 



where 

2 o(u) _ 2 

— j Pi — ; 
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pt } = Z> Pi U) = -^- (6-9) 



Comparing Eqs. (P9|) and ( |S.9[) one sees at once that the beta-functions coincide up 
to two loops if one could identify u with g. On the basis of the naive continuum limit 
( |6.4|) this is of course not possible as it would rather suggest u oc g/2 corresponding to 
the (nonchiral) GN model. The perturbative free energy Eq. ( pid| ) contains, however, 
the two loop beta-function of the CGN model already as an input, and since Eq. (P10| ) 
is in perfect agreement with the result obtained in the scaling limit of the Hubbard 
model ( |5.9| ) this provides a rather strong (though indirect) evidence that the two loop 
beta-functions do in fact agree. It would be clearly desirable, however, to show in a 
more direct way the equivalence of the beta-functions and understand to some degree 
the problems encountered after having taken the naive limit ( |6.4j ). 

The one-loop renormalization group analysis outlined in Appendix |F| gives an indication 
that the rather pathological theory obtained in the naive limit (|6.4j ) would eventually flow 
to the CGN model. (We note that this result has already been indicated by the "g-ology" 
treatment of the ID electron gas |35fl .) As the renorm trajectory of (|6.4|) approaches that 
of the CGN model the couplings become large and the one loop analysis breaks down and 
that is why the one loop analysis is not conclusive. 
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APPENDIX A: 



In this Appendix we would like illustrate our proposition about the isospin structure of 
the limiting states on an explicit example. 

Suppose we have a certain (77, i?) set with n r = 3 solving Eqs. (|3.9| a-b) and a (£, <p) set 
with rii = 2 satisfying Eqs. ( |3.1(J| a-b). In the following we enumerate the possible isospin- 
structures of states characterized by such a set (77,$) and (£,</?)■ F° r the considered values 
of n r and ni there are three possibilities: there are no 9s, there is one 9 (which has the value 
9 = -2/5 due to (|3TT3|)), and there are two 9s (whith values 9 ± = -2(1 ± i^2)/3). This 
way the the BA eigenstates (longest projection states) are 



(Al) 



(Here the two numbers refer to the length and the third component of the isospin.) As we 
argued, based on the independence of the 'r' and T sectors it is natural to think that if 
n(9) = both sectors are longest projection isospin eigenstates and the state factorizes: 



»7> f, 


5\ 
2/ ' 


n(9) 


= 




3 \ 
2/ ' 


n{9) 


= 1 


»7, 5, 


1 \ 

2/ ' 


n(9) 


= 2 



77, 0;£,p;§, § 



(A2) 



(Although this is a natural assumption, it is not trivial: before taking the scaling limit the 
two sectors are not separated, so they and their isospins are not defined; while in the scaling 
limit the wavefunction does not exists, so ( |A"2j ) can not be directly checked. Nevertheless the 
assumption, that the isospins of the two sectors in the scaling limit can be defined separately 
leads to a consistent picture.) Applying I~ on ( |A2| ) we get a state 



3. 3 1 

' 2' 2 



According to our proposition the state with one 9 is of the structure 

77, 7? ; e, 1 1) = B lfi \ v , |, |)je, ^ 1, i), + B 0tl \n, #] §, |)Je, i, o 

Acting on ( |A3a| ) and ( |A3b| ) by I~ leads to 



(A3a) 



(A3b) 



77, 79; 5 



2.0 



' 2' 2, 

^ 0; i> -i)J^> !> i), + a ia\v, #i §, w 1, o) { + A 0)2 ^ tf; |, §)Je, <^; 1, -1) 



(A4a) 
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and 



' 2' 2 



B 



2,0 



V, I -|)Je, v>; l, l> ( H- Bi,i\v, f , |)Je, l, o) z + S ,2^ f , f l, -i),. 

(A4b) 

At the same time the state with two 6s is of the form 



a 



2.0 



77, •& 



' 2' 2 

■ 3 _i 

' 2' 2 



£, 1, 1), + Ci,! 77, 7?; §, £, 1, 0) j + C , 2 77, 0; f , |) r £, 1, 

(A4c) 

By further applications of J - we get 



.4 



3,0 



' 2' 2 

v, I -l\ f , !> i), + ^ ^ ^ ; I' -0, ^ ; x ' °); + Al > 2 ^ ^ ; I' 0. £> ^ _1 >,' 



(A5a) 



5 3 ,o »7, 0; |, -1) j£, ¥>; 1, l) + B 2A \ v , 0; §, -|)Je, 1, o) j + £ 1)2 |t7, 



■ 3 I 

' 2' 2 



(A5b) 



a 



3.0 



77, 7?; |, -|) r £, 1, l) f + C 2)1 77, 7?; |, -1)^ £, 1, 0) ; + C lj2 77, 0; |, ^ ? 



(A5c) 



77, 7?; £, f , -I) = A 2j2 |t7, 7?; §, -§}j£, ^ 1, + A 3 ,i|r/, 0; |, -§)j£, <p; 1, o) f , (A6a) 

^ ^ 1> -|)J^ V5 !> + fi 3,i|r/, 0; §, -f)j£, 1, 0^ ; (A6b) 



,.3 3 



v, 2 ' 2 

and 



(A7) 



(The various coefficients should be proportional to the corresponding Clebsh-Gordan coef- 
ficients.) All these states are of the same energy and momentum. In the space spanned by 
them we can choose for a basis the states 



77,7?;|,J r 3 )Je,^l,/^ It 



■>\ r> - § I _i / 3 = 1 n _i 

2'2' 2' 2' 1 ' ' 



(A8) 



These basis vectors can be expressed in principle through eqs. (|A2|)-( |A"8|) . It is obvious, 



that this procedure can be performed for any (77, £, tp) solution of Eqs. ( 3.9a , 3T9rj ) and 
( p.l0a| , p7lUb| ), thus the space spanned by the scaling limits of the solutions for the BA 
equations factorises into the direct product of the spaces spanned by the isospin eigenstates 
of the 'r' and T sectors. 
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APPENDIX B: 



Deriving the HLBA equations one neglects corrections exponentially small in the chain- 
length (terms of the type 5 oc exp{— a(u)N}) at two points: 1) when describing the bound 
pairs of ( [2.4|) by a single parameter (the A), and 2) when representing the complex As by 
the set 0. Although in the scaling limit — > oo, a(u) may tend to zero, thus one has to 
check if these corrections remain negligible indeed in the scaling limit. This is what we do 
in this Appendix. 

1) First we discuss the correction terms to the size (in wavenumber space) of the bound 
pairs in the scaling limit. The relation fl2.4p completed with the correction terms reads 



sink* = A± iu + 5 ± . (Bl) 



In (|BT| ) 5 can be estimated in the following way. Substituting (|B1|) into (|2.2a|) we have (up 



to leading order in 5) 

?l I A -i- in — A' 1 

(B2) 



In — = -Im { Nk + - V 2 tan" 1 } 

2U { A'(#A) U J 



The sum over the A' (in leading order in N) can be calculated by means of the ground-state 
density of A: 

l n ffi = _ivim \k + - J 2 tan' 1 A + m ~ A ' a JA') d A> 1 ( B 3) 
1u [ J u J 

where (according to [§]): 

a oW = ±] J °M™^ <L,, (B4) 

27T J cosh I CUM I 

v ; 

with Jq being the zeroth order Bessel function. As for « < 1 the main contribution of the 
integral comes from the region |A'| < 1, evaluating the r.h.s. of (|B3| ) the regions |A| > 1 and 
|A| < 1 have to be treated differently (for the sake of simplicity we consider here real As 
only; a complex A can be treated in a similar way) . In the first case for u 1 the integral is 
real, (sgn(A)7r/2) giving no contribution to the imaginary part, and then the task is simply 
to calculate the Im k + 

Im k + = cosh" 1 1 - (^u 2 + (A + l) 2 + v /W(A-l) 2 )j . (B5) 
For the second case the integral can not be neglected, and we have to calculate 
cosh" 1 ji (\/W(A + l) 2 + \J v? + (A - l) 2 ) | - 



1_ J ? (A - A') 2 JpMexp^AQ 
8tt J J Au 2 + (A - A') 2 cosh(cuu) ' 1 ' 



— oo — oo 
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This expression can be evaluated (approximately) for small u. Finally we obtain 



m f-JVcosh-^AI (|A|>1), 

with Kq being the modified Bessel function. This shows, that in the scaling limit ( |1.2|) 5^0 
everywhere except in a small vicinity of A = 0. There it behaves as 

\5\ = 2-uexp j— LniQ cosh | . (B8) 



Consequently the bound pairs are of the form ( |2.4| ) indeed up to exponentially small cor- 
rections. Moreover, at the ends of the A distribution, where the holes corresponding to the 
massless particles are situated, Q2.4p becomes exact in the scaling limit. 



It is apparent that in the scaling limit the 5s in ( |B8|) can be larger than 1/N. Nevertheless 
their effect is negligible compared to the leading finite size corrections, which are actually of 
the order of 1/L (see Sec. [TV] and Appendix |Cj), as the coefficients of the 1/N terms are of 
the order of N/L. 

2) Now we discuss the neglected corrections to the complex As. The complex As are 
represented by the set in the following way: a with an imaginary part |ImO| > u 
represents a wide root: 

— ► A = + msgn(Im0) (|Im0| > u) , (B9a) 

and a with an imaginary part |Im0| < u represents a close pair: 

— >A ± = Q±iu±r (|Im0|<u), (B9b) 

with r expected to be exponentially small. Its magnitude can be estimated on the basis, 
that A 1 * 1 must satisfy Eq. (|2.6p. Substituting A + into (|2.6|) we have 



2nJ = N (sin~ 1 (A+ - iu) + sin" 1 (A+ + iu)) - ^ 2 tan" 1 - A WU * 



k u 



,A+-A' , A+ - A" , 

- Y 2 tan" 1 2 tan" 1 . BIO 

^ 2v 2?/ 

Neglecting the effect of the ks (as being not macroscopic) and estimating the sum over the 
A' by means of the ground-state A distribution a , ( B10| ) yields to leading order in r 

\t\ ( / \ r°° A + A' 1 

- In Li = Nlm \ (sm-\A + - iu) + wrHA + + iu)) - 2 tan" 1 a (A) \ . 

2u v ' J-oo 2u 



(Bll) 



The r.h.s. for small u can be evaluated in a closed form if the |ReA + | > 1. The result for 
ReA + > 1 reads: 

7rImA + / 7r \ f 7rReA+ ] 

iv2sm — j »y exp {-— }■ < bi2) 
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which in the scaling limit yields 



jp = exp {-Lm cos had e Re ®} , (B13) 
with $ (in agreement with (|3.7|) ) 

*=-« A+ - iu - i K (B14) 
The analogous expression for the ReA + < — 1 (which are connected to the T particles) 

- — - = exp j— Lmo coslm^ e~ Rcv j , (B15) 

with 

7r(A+ - m + 2) . , 

^=-^-2= (B16) 

The conclusion of the above is, that the corrections to the size of the close pairs are 
negligible indeed as long as Re$ 3> — In (Lmo) an d Re<£> <C ln(Lm ). Ret? ~ — ln(Lm ) 
and Rey? ~ ln(Lmo) can occur only in the smallest energy states (when the r] ~ — ln(Lmo) 
and £ ~ (Lmo)), i-e. one has to treat only the smallest energy states in a more subtle way, 
all the others are correctly described by Eqs. ( pT9| ) and ( |3.10| ). 



APPENDIX C: 

In this Appendix we would like to discuss the finite size corrections in more details. The 
procedure given in Sec. ^ is based on the fact, that the energy-momentum dispersion rela- 
tion of the massless (dressed) particles is linear, and that the dressed particles are described 
by BA type equations. Here we describe the better known method [[JI],[8|, which does not 
use these properties directly. 



Our starting point is given by the Eqs. ( p.5| )-( [2l)|) and the dispersion relation (|2.7| ). For 
the sake of simplicity we consider states with no massive particles. The A distribution 
is supposed to be the simplest, i.e. the J v parameters form an equidistant series between 
J + <(N- n(A) - l)/2 and J- > —(N - n(A) - l)/2 (with J+ - J- = n (A) - 1), i.e.: 

J = Ji, J v + l — J v+ x, J n (A) = J + ■ (CI) 

After introducing the function 

if/ \ 1 n(A) A - A 1 

2(A) = — (sin-^A - iu) + sin-^A + in)) --^ 2 tan" 1 , (C2) 

through the equation 

JJN = z(A v ) (C3) 
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one can interpret the A,, as a function of J v . The density of the roots is given by 



Defining 



■KA) = -$>. (04, 



S{A) = ^6(A-A„)-*{A), (C5) 



we can write 



If//, • ^ 2 \-V2 / 1 , . v 2 \-i/2 



*(A) = — l((l-(A-iii)') +(l-(A + i« 



'-oo 4u 2 + (A-A') 
This equation, after some straightforward algebra leads to 



4it 1 
' -HA') + S(A')) . (C6) 



a(A) = (T (A) - J K(A - A')S(A') , (C7) 

— oo 

with cr given by ( |B4| ) and 

R(A) = — / „, , . (C8) 

— oo 

The energy 

E = -Nu - 2t (y/l - (A v - iu) 2 + yjl - (A v + iu) 2 - 2u\ , (C9) 
jj \ / 

after some manipulation turns out to be 

+oo 

E — Ne — J e c (A)S(A) , (CIO) 



with e an d £ c given by Q2.ll ). To calculate integrals of S(A) the Euler-Maclaurin type 
formula 



I J 1 n_1 J I J Jn ^ N 1 J, 

2N f \w) + N J/ ( Nj + 2N f ( if) = I f ^ + 12l^V\W)~ f 



Ji/N 



(Gil) 



(/' being the derivative of /) can be used. (We note here, that ( |C11|) applied in this 
and related problems (like the Heisenberg chain) takes into account the leading terms of 
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a nonconvergent (asymptotic) expansion. Nevertheless the result is reliable, as it is well 

(C12) 



indicated by the numerical calculations for very long chains [30].) This leads to 

E — Ne = Ne + + Ne~ , 
where 



A+ 



e(AMA) - ^(A+) - ^f^L 



(C13) 



with A + being the largest root: 



A+ 



(C14) 



and analogous expressions give e~ and A - . It is apparent, that to calculate e + , it is enough 
to know cr(A) for A + < A < +00. Applying (|C11| ) in (|C7|) and neglecting the next to leading 
terms one arrives at the equation 



o-(A) = o-o(A) + / K(A - A'W(A') - —K(A - A+) + K '^\ K+ ] 



(C15) 



wich determines <r for A > A + with sufficient accuracy. ( |C15|) is a Wiener-Hopf type 
equation, and can be solved in closed form |30|j38|1 leading to 



£+ = 8vrt 



1 (n 



+^2 



The analogous expression for e~ is 

M%) J 1 



8 N 2 48N 2 



n 



A" - n(A) - 1 



(C16) 



8vrf 



n 



1 



A" 2 48N 2 



n 



N-n(A) - 1 



J - 



(C17) 



The momentum P = J^/N through an elementary calculation is found to be 

o 77 / + 2tt f(n-) 2 (n+) 2 l 

P = - n ) H <^ - — i - — i- ^ . 

2^ AM 4 4 1 



(C18) 



Completing the scaling limit, and renaming n one arrives at 



E-Ne = -— + — 
bL 



7r 27T J (n r ) 2 (n^) 



+ 



2ir \ (n r 



4 



(C19) 



This is the same as obtained in Sec. [IV]. 
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APPENDIX D: 



Here we outline the main results of the perturbative computation of the free energy in 
the CGN model p3]| . First one introduces a chemical potential, h, coupled to a Noether 



charge of the U(2) symmetry, Q, and considers the ground state of the system with the 
Hamiltonian 

H = H -hQ. (Dl) 

The charge matrix Q has the eigenvalues Q a on the bare fields: Qip a = Q a ipa- By decom- 
posing this diagonal matrix into abelian (SU(2) singlet) and non-abelian (traceless) parts 
one can write: 

Q ab = qS ab + Qrl b T*= l -a\ (D2) 

with a % (i =1,2,3) being the usual Pauli matrices. As mentioned above, in the CGN model 
the massive particles couple only to the non-abelian part of the charge, Q, while the massless 
sector couples only to the abelian component, q. Consequently, the free energy (the ground 
state energy of (pi|) ) is expected to have the form: 

8f(h) = f(h) - /(0) = f(h, q, g') + f(h, Q, g(h)) , (D3) 

where g' = g v +g s /2 is the coupling of the abelian sector. Perturbation theory is applicable for 
weak coupling, which implies high densities (asymptotic freedom) as the running coupling, 
g(h), behaves as g(h) cx [ln(/i/A)] _1 , where A is a parameter (of dimension mass) which is a 
renormalization group invariant combination of the ultraviolet cutoff and the bare coupling. 
The free energy of the massive sector, /, has in fact the form / = h 2 Q(h/A,Q) and by 
expressing the running coupling, g(h), in terms of the A-parameter one obtains an asymptotic 
series in [ln(/i/A)] _1 for $(/t/A, Q). As the Abelian coupling g' has a vanishing beta function 
f{h,q,g') is simply given as h 2 (p(q,g'). The calculation of the free energy in the scaling 
limit of the Hubbard model based on the HLBA eqs. yields a result of the form 5f(h) = 
h 2 ty(h/m ) + /i 2 const, and when for a large density f(/i/mo) is expanded into an asymptotic 
series in l/[ln(/i/m )], the results for $ and ^ should agree as one calculates the same 
quantity provided the two models are indeed equivalent. One needs however to relate the 
two mass scales mo and A in order to achieve agreement. 

The abelian part of the free energy, /, receives contribution only from bubble diagrams 
whose sum yields: 

h 2 a 2 1 

/ = "— T-V7T- (D4) 
71 1 + 2g'/7T 

This can be immediately seen to agree with the result ( |5.14j ) found in Section V for the 
massless sector by identifying the chemical potentials provided that we set g' = 0. 

The non-abelian part of the perturbative free energy had to be calculated up to three 
loop order with the result: 



1 2tt 2 



9(h) . „9 2 (h) 



+ C^ + 0(g b ) 



IT 7T 2 



(D5) 
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where g(h) is the renormalized running coupling in the MS scheme and the constant C is 
given as 

C = 21n2-- + 21nQ. (D6) 

Introducing now the renormalization group invariant Aj^g— parameter corresponding to the 
running coupling: 

(h) = 1 _ A ln\n(h/A m ) / \n\n(h/A m ) \ 

9{ > (3oHh/A m ) ^(MVAms)) 2 V(MVAms)) 3 / ' 1 ' 

where /3q and f3± are the universal 1- and 2-loop beta-function coefficients. The beta-function 



of the SU(2) CGN model up to two loops is given as |37 

k §h 9{h) = ~^ g2 ^ ~ " ' • • ' ( D8 ) 

A>=-, Pi = -- 2 - (D9) 

In order to facilitate the comparision of the perturbative free energy with the nonperturbative 
result derived from the HLBA eqs. of the Hubbard model we rewrite the perturtbative free 
energy flB7|) using \n(h/A m ) = t + \nR m , where R m = m /A m : 

JK J 27r 2 I 2t At 2 It 2 I t 3 / I K ' 



where 



C=-C + \nR m (Dll) 



APPENDIX E: 

We sketch here briefly the 'naive' continuum limit of the half filled Hubbard chain. A 
more detailed description of the procedure is given in Ref. [BU . 



Let us redefine the chain so, that four lattice sites form one elementary cell and define 
the operators 

1 4 

</>vM = TrT,(- i Y 3c M»-i)+w " = 1,2,3,4, a =U . (El) 

4 V a j=i 

In terms of these the original Fermion operators are 

4 

CA(n-i)+i,o = V^EW^^W Z = l,2,3,4. (E2) 
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The length of the chain is 

L = Na = AaN' (N' = int.) 
and the continuum limit is defined as 

a — > 0, N' — > oo, L = fixed 

with the continuous variable 

x = 4a(n - 1/2), dx = 4a. 

In this limit 

^ f dx 
^ ^ 7 4a 



(E3) 



(E4) 



(E5) 



If we introduce 



8n,ri 4a5(x — x). 



0i, CT (n) = i/)i,ct{x), h,a(n) = iJ 2 ,a(x), 



(E6) 



(E7) 



{^(x),^,^')} = 5(x - x')8 ati3 8 ai(T >. 



(E8) 



Applying (|E2]), ( |E7| ) and ( |E5| ) to the Hubbard Hamiltonian and neglecting 2 and </>4 
(this can be done as these fields become infinitely massive in the continum limit) one obtains 
a Hamiltonian which is not even SU(2) symmetric. Observing, however, that the interaction 
part of the Hamiltonian ( | . 1|) can also be written as 



H, 



U 



N 



N 



int 



H'-UN + -UN where H' = — J2( n i1 + n jif and N = Y.( n jt + n jl) 



one obtains in the naive continuum limit the following Hamiltonian density: 
H{x) = (2a*) (- £ ^75^ CT + u [(E^7o^) 2 + 

where 



^(7 



1, (7 

2, (7 



75 



i 
-i 



7o 



1 

1 



(E9) 



(E10) 



(Ell) 



and u = U /At. (Note that this differs from the definition used in the bulk of the paper as 
there u = \U\/±t.) 
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APPENDIX F: 



In this last Appendix we would like to present a one loop renormalization group analysis 
of the naive continuum limit of the half filled Hubbard model fl6.4p . 

Let us consider the following non Lorentz invariant theory with four fermion interactions: 



£ = iip a 0ip a 



9o, T 



910/ T 



#11/ T 



92, J 



(Fl) 



The theory given by ( |bTp is clearly renormalizable at least up to one loop and contains 
the theory obtained in the naive continuum limit ( |6.4| ). The CGN model now corresponds 
to the choice of 0o = 02 = ^9io — 2gn, while the nonchiral GN model is obtained as 
92 — 9w — 9n — 0. The theory found in the naive continuum limit (|6.4|) corresponds to 



9o = 9io = 2n, g 2 = gu = 



(F2) 



Our strategy is to compute the one loop /3-functions in the space of the four couplings 
{<7o j 0io j 5*11 j 92} an d then investigate if the classically neither Lorentz nor chiral invariant 
theory ( |6.4|) could somehow be connected to the SU(2) CGN model in the one loop ap- 
proximation. The calculation of the one loop /3-functions in perturbation theory is by now 
standard and therefore we just state the result of the calculation without going into the 
details. 



0o 

02 
010 
011 



1 


(n — 


7T 


1 


(n — 


7T 


1 




0002 , 

7T 


1 




9o92 , 

7T 



n - l )gl + 9o92 - (910 + 9n)(9o ~ 92, 
n ~ l )gl + 9o92 + (910 + 9ii)(9o - 92) 



(F3) 



where 0i = [idgi/dfi. The calculation has been performed for the SU(n) symmetric case 
in the dimensional regularisation scheme. One avoids the problem of evanescent operators 
since in the computation of the /^-functions up to one loop order, only the divergent pieces 
are needed hence one can perform the Dirac algebra in two dimensions. For the Lorentz 
invariant case gio = gu Eqs. ( |F3| ) reproduce the known result of Ref . JHJ . For the present 
case of interest n = 2 and then one can easily solve Eqs. ([F3| ) for arbitrary initial conditions. 
Starting with the initial conditions ([F2j) at t — the evolution of the couplings is given as: 



9o = u 
92 = u 

u 



1 1 

+ 



l + 2ut/n l-2ut/n / 

: 1 1 ^ 

l + 2ut/n l-2ut/ic t 
( 1 1 



(F4) 



010 

0ii 



2^1 + 2ut/n 1 - 2ut/n / 

u ( 1 1 
+ 



2 V 1 + 2ut/n 1 - 2ut/n 



u 



u 
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Independently on the sign of the coupling of the Hubbard model, u, the couplings 
{<?o> 9i0i 9iii 92} w iH inevitably become large for t — > ±00 (both in the ultraviolet (UV) 
and in the infrared (IR)) and perturbation theory then breaks down. In other words one 
hits a Landau pole both in the UV and in the IF region which appears to be rather peculiar. 
For u > in the t — > 00 limit the couplings will approximatively satisfy 

goSiusi-^Ri 2g 10 « 2#ii , (F5) 

while for u < 



g « u « 5-2 ~ 20io ~ 2#n . (F6) 



In the IF limit the situation is reversed in that then one gets close to the CGN trajectory 
for u > 0. Unfortunately the one loop beta-function of the Lorentz-invariance violating 
interaction (gio — gu) is zero, so in the one loop approximation no definite statement about 
the restoration of Lorentz invariance can be made. Clearly it is tempting to conclude from 



Eqs. ( p5| , p6|) that Lorentz invariance will be restored and that the IF limit of theory (p.4|) 
for u > resp. its UV limit for u < will eventually be the CGN model. As the one loop 
approximation is applicable, however, only for t l/l u l> the above can only be considered 
as an indication about the limit of ( |6.4|) . Nevertheless it is quite interesting that this might 
be an explicit example of the restoration of Lorentz invariance in the quantum theory. 
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